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Abstract 

The dynamic transition between the ordered flow and the plastic flow is 
studied for a two-dimensional driven vortex lattice, in the presence of sharp 
and dense pinning centers, from numerical simulations. For this system, which 
does not show smectic ordering, the lattice exhibits a flrst order transition 
from a crystal to a liquid, shortly followed by the dynamical transition to 
the plastic flow. The resistivity provides a critical order parameter for the 
latter, and critical exponents are determined in analogy with a percolation 
transition. At the boundary between a pinned region and an unpinned one, 
an anomalous layer is observed, where the vortices are more strongly pinned 
than in the bulk. 
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I. INTRODUCTION 



Following extensive studies on the effect of disorder on the static vortex lattice, the 
physics of the vortex lattice with random quenched disorder and driven by a uniform force 
has attracted recently much attention. Interacting systems, forming periodic structures at 
the equilibrium, were already the subject of much interest since the earlier studies of charge 
density waves The complexity of the depinning phenomena was soon pointed out, 

in the sense that the description of the depinning threshold by a critical phenomena is no 
longer valid when one takes into account the possibility of topological defects within the 
periodic structure Plasticity, which is commonly observed closed to the depinning 

threshold, is a dramatic illustration in the case of the two dimensional vortex lattice. There 
has been several investigations of the driven vortex 'phase diagram' which have enriched 
the canonical description [Q: pinned vortex glass - plastic flow - moving crystal as the 
driving force is increased. Amongst these, numerical simulations of two dimensional vortex 
assemblies, initiated by the work of Brandt 0], have very often accompanied theoretical 
progress on the subject. After the proliferation of the lattice defects was put into evidence, 
suggesting a dynamic first order melting transition at the occurrence of plastic flow 0,0, 
numerical simulations identified the ordered phase as a moving transverse Bragg glass [0, 
in agreement with theoretical expectations pi. Simultaneously, Ryu et al showed that an 



hexatic order parameter exhibits a sharp transition as one enters the moving glass [|lT 



Later, Olson et al showed that, for a soft flux lattice {Ay < 1), the ordered phase presents a 
smectic order |12|. Finally, Kolton et al recently introduced a 'frozen transverse solid' beyond 
the smectic regime, characterized by a drop of the Hall noise [0. In a general way, there 
is often some confusion about the exact nature of the 'transition'. 'Dynamical transition' 
and 'phase' are often employed in place of 'crossover' or 'dynamic regime', without further 
justification. Indeed, there seems to be up to now only one strong suggestion of a genuine 
dynamical transition in the works in refs [0] and [|ri[] . The fact that the notion of dynamical 
transition is itself defined only with difficulty (see ref. p|) has certainly contributed to this 
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situation. It is not clear, for instance, wether one should try to use some dynamical quantity 
- such as the correlation length of the local velocities - as was done in ref . , or if one should 
use some instantaneous, topological one - such as the concentration of defects or the hexatic 
order parameter in refs |Q and - in the search for an order parameter. Here, it is shown 



that a simple system, not showing any intermediate smectic order between the ordered and 
the plastic flow, exhibits a second order like transition to the plastic regime. 

II. EXPERIMENTAL DETAILS 

A two dimensional lattice subjected to a uniform driving force (applied along the y-axis, 
thereafter denoted longitudinal direction) in the presence of pins is simulated, using the 
force equation : 

fvv(r) + fp(r) + fBo(x) + J A $0 - r = (1) 

The geometry is analogous to the one of a Corbino disk experiment: the two edges at 
y = const, are submitted to a periodic boundary condition; the ones at x = const, are 
submitted to an external magnetic field, Bq, which is simulated by an extra force fsg acting 
on each vortex, perpendicular to the edges. The forcejBoix), acting on a vortex at a distance 
X from the edge, is that imposed by a semi-infinite vortex lattice at a distance + x, where 

1 /2 

ao = ($o/-S) is the flux lattice spacing at the equilibrium. Flux lines are assumed rigid 



rods and the force per unit length between vortices separated by a distance r is [jT3[: 



fUr) = {Ay/\)K^{r/\) (2) 

where Ki is a Bessel function. This is strictly a good approximation only in the case of vortex 
lines (rods) and for 2D vortices a logarithmic interaction should be used. The interaction 
between vortices was cut at a distance 5A, using an interpolation to zero. This was done in 
order to avoid spurious distortions of the equilibrium lattice from the Abrikosov lattice or 



the introduction of topological defects, as was shown to occur for a sharp cutoff in ref. [IT3 . 
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The sample dimensions were 100 along x-axis and 70 oq along y-axis. Strong pinning 
centers are randomly distributed in the sample. A pin free region was left for x < 25 ag and 
X > 75 ao- Doing so, a defect free lattice is obtained at the edges of the sample, providing 
well defined boundary conditions. The density of the pinning sites is np = i?#/$0; with 
$0 the flux quantum and the 'matching field' for which an equilibrium flux line lattice 
shows the density of flux lines ny- The force per unit length exerted by a pin at a distance 
r from the line is given by: 

/p(r) = (2 Ap/rp) (r/rp); for r < rp, for r > rp (3) 

The pinning force is exactly balanced by the Lorentz force, 3 A^o, for J = Jq = 2 Ay/rp $o 
(in the following, j = J/Jq). In the present study, the following parameters were used : 
A = 1.57 ao, rp = 4.9 10"^ oq, Ap/Ay = 2.5 10"^ and B^ = Q B. Using the notations in 0, 
this corresponds to sharp, dense and strong {Ap/r^ oq cee ^ 1) pins. The sample contained 
approximately Ny = 7000 vortices and 25 000 pins. 

III. RESULTS AND DISCUSSION 

Vortices trajectories are shown in Fig.|l]. At first sight, they display a striking feature: 
as the driving current decreases and the trajectories evolve from correlated channels to 
branched trajectories, the vortices are first pinned at the interfaces between the pinned and 
the unpinned region. This is in contradiction with the intuition gained from fluid dynamics 
physics, where one would expect the average velocity of the fluid to decrease monotonously 
from the one for unpinned vortices to the one of vortices slowed down by solid friction. 
Rather, as shown in Fig.^], the average velocity first drops to a minimum right at the interface 
between the pinned and the unpinned region, and then grows to some roughly uniform value 
at a distance ~ 5 ao from the interface. The magnitude of this anomalous boundary layer 
effect may be measured as the ratio of the average velocity in the layer, to the one far 
away in the pinned region (Fig.^c). Dynamics regimes were characterized using physical 
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quantities as commonly done in flux lattice simulations [l7|,p^p|. As shown in FigJ^, the 
system exhibits a sharp departure from a linear V — I characteristic; an onset of the voltage 
noise measured in the direction transverse to the average flux flow; an onset of the lattice 
diffraction peaks widening as well as the onset of the anomalous layer effect at ji 0.33. 
Close to this value, at j2 — 0.305, the voltage derivative, dV/dJ, shows a sharp peak; 
diffraction peaks vanish and the layer effect saturates. The analysis of the structure factor 
S(k) = \ e'**^""'!^ on an annulus which overlaps the first Brillouin zone diffraction 
peaks (Fig.^ shows the progressive evolution of the central region of the sample from a well 
ordered hexagonal lattice at j = ji to a liquid at j2- In between, there is no evidence in 
the diffraction intensity for an asymmetry between the average flux flow direction and the 



one transverse to it. Following ref. [12|, the regime at j < j2 is that of the plastic flow of 



the amorphous solid. The transition region j2 < j < ji between the plastic regime and the 



ordered state differs from the ones described in |jT2| or 0, as we find no evidence for the 
asymmetry needed in the diffraction intensity for a smectic order or an order intermediate 
between a smectic and a crystal. Also, the transition regime width observed here is only 
about 10% of the critical value for the plastic to quasi-ordered regime current, while values 
larger than 30% were found in |T^. These differences are due to parameters much different 



from the ones used in [|l2l. Here, pinning sites are dense and almost point like {np/ny = 6 



and ab/rp = 25), while the pinning density is comparable to the vortex density and pinning 



sites are extended in |jT2[ (np/n^ = 1.4 and ab/rp = A). As a result, vortices do not sense 
here the asymmetry of the pinning potential (when it is tilted by the driving force) as they 
do for extended defects, and the smectic regime does not occur. 

Within the plastic regime, the evolution of the channels resembles that of a percolation 
transition, and the transition between the ordered flow and the plastic regime may be viewed 
as the percolation of dynamic flux channels in the transverse direction. This similarity was 
already noticed earlier in |]1T|. The analysis of the resistivity quantitatively demonstrates 
the validity of a critical phenomenon approach. As seen in Fig.^fe, the resistivity may be 
fitted to a critical order parameter of the form (1 — j/jc)^, with (3 = 0.34 ± 0.02. The 



critical driving force obtained in this way, jc, is within fitting uncertainty identical to j2- 
The restricted intermediate regime, as observed here, might be crucial for the observation 
of the critical behavior, as it tends to smear out the transition. The interpretation of the 
anomalous layer effect - which is fully developed once one has entered the plastic dynamical 
phase (as defined from the critical analysis above) - appeals for a better understanding of 
the latter. Characterization of the instantaneous structure, such as the structure factor 
displayed above, is useless to the study of the plastic phase: the autocorrelation function of 
the instantaneous vortices positions , C(k) =< p(r) p(r + k) >r, where p(r) = 
only confirms an evidence for a liquid order (Fig.^. Considering the existence of two distinct 
vortices populations [0 : a rapidly moving ensemble of vortices along quasi static channels 
and quasi pinned ones, one may also define the velocity- weighted autocorrelation function, 
Cv(k) =< pv(i") Pv(i" + k) >r where py(r) = (^(i^i) This essentially measures the 

correlation amongst the most mobile vortices. As can be seen in Fig.^ the function evolves 
from the one of a liquid to the one characteristic of isolated flux channels (two peaks at 
k = (±ao, 0)) as j decreases. Both C(k) and Cv(k) show that second neighbor correlations 
are strongly damped in the transverse direction. However, correlations between vortices 
may be found that are less demanding than the ones uncovered by the transformations 
of the instantaneous lattice. The autocorrelation function of the channels, defined as : 
Cc(k) =< pt(i') Pt(r + k) >r where Pf(r) = /Qpv(r) dt and t is a time large enough so 
that the most mobile vortices have moved by a distance larger than oq, provides evidence - 
close to the transition - for stronger transverse correlations between such channels than the 
ones uncovered by C(k) or Cv(k) (Fig.^. This means that channels tend to correlate in the 
direction transverse to the average flux flow. 

The qualitative analogy with percolation and the definition of a critical order parameter, 
as shown above, both appeal for a definition of dynamic clusters. This is done in the 
following way: first, pt(r) is computed as defined above, thus providing some snapshot of 
the channels. Then, the pattern defined in this way is filtered from frequencies larger than 
Oq ^ (this insures that two contiguous channels do belong to the same cluster). Finally, one- 
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dimensional clusters are defined, as the line segments perpendicular to the average flux flow 
that are entirely contained within the filtered channels (Fig.^. Such a definition takes into 
account the anisotropy of the problem and insures that an infinite cluster is found at the 
transition. Although it provides an infinite number of clusters for each sample, it allows the 
study of clusters distributions as commonly done in the study of percolation [0. The first 
infinite cluster is found at j = j2, in agreement with the critical analysis of the resistivity. 
As shown in Fig.^ it is found that the mean cluster size, S = J^s^"^ scales close to j2 as 
(1— j)"''' with 7 = 1.2±0.02. For j < 0.2, the mean cluster size saturates to S* = a^, meaning 
that one enters a regime of isolated flux channels. This agrees with Cc(k) in Fig.^, where it 
is seen that the first neighbor correlation roughly become isotropic below j ~ 0.2. Following 
the analogy with percolation, one may also define a correlation length, ^, which diverges 
at j2- Then, the saturation observed for S may be directly interpreted as the decrease of 
^ down to the average flux line spacing, oq. Within this description, it could be appealing 
to interpret the existence of the anomalous interface layer as a 'proximity effect'. However, 
the order parameter - as defined above - should in this case continuously increase from zero 
in the ordered phase, to the value of the bulk over a distance comparable to ^, whereas it 
is anomalously large in the pinned layer. Also, the thickness of the anomalous layer should 
strongly depend upon j, which is not observed in Fig.^ A more plausible interpretation for 
the effect is in fact a topological one. Channels transverse wandering - an alternative view for 
the clusters distribution and the fractal topology of the plastic phase - is strongly suppressed 
at the interface with the ordered phase. Besides the occurrence of the ordered phase, the 
occurrence of a pinned region at the interface provides another way to pin the transverse 
excursions of the vortices as imposed by the proximity of the crystal structure, hence the 
observed effect. However, this piece of explanation does not provide any estimation for the 
width of the layer. This shows that, although the present analysis provides some evidence 
for the existence of a dynamical plastic phase and a second order like transition, we still lack 
a complete understanding for the pinned, driven vortex lattice. 

I gratefully acknowledge valuable help from S. Ravy in the handling of numerical diffrac- 
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FIGURES 





FIG. 2. Average y-velocity component profile. The pinning centers density is non zero where 
a; > 0. For clarity, results for each curve were rescaled along the vertical axis. 
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FIG. 3. a, filled : average voltage along the main flow, computed in the pinned region, away 
from the anomalous layer; a, empty : voltage derivative. 6, points: difference between voltage 
in a and the free flux flow voltage; b, line : fit to Vb(l — j/jc)^, where jc = 0.304 ± 0.001 and 
P = 0.34 ± 0.02. c, crosses: transverse Hall noise; c, circles: average velocity in the anomalous 
layer, normalized to that of the bulk; c,triangles: inverse of the width of the diffraction peak at 
k = (ao, 0). The vertical line at j = 0.304 marks the critical driving current, as given by the fit in 
6; the one at j = 0.33, the onset of departure from the free flux flow potential in a. 
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FIG. 4. Polar plot (logarithmic units) of the instantaneous structure factor, after radial inte- 
gration over an annulus which overlaps the diffraction peaks in the first Brillouin zone (the flux 
lattice is sampled in the pinned region, away from the anomalous layer). The zero angle axis points 
along the reciprocal direction transverse to the average flux flow. 




FIG. 5. Gray scale maps of autocorrelation functions. Clockwise: j = 0.31, 0.295, 0.27, 0.24. 
Left upper quadrant: autocorrelation function of the instantaneous lattice, C(k). Left lower quad- 
rant: autocorrelation function of the instantaneous lattice, weighted by the y- velocity component, 
Cv(k). Right half : autocorrelation function of the vortices trajectories, weighted by the y- velocity 
component, Cc(k) . The average flux flow is horizontal. 
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FIG. 6. Channels structure {j = 0.27), filtered from frequencies higher than Qq . Shown as a 
white line is a one-dimensional cluster of size s. 
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